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Abstract. Ultrafast optical pump-probe spectroscopy measurement on monolayer 
graphene observes significant optical nonlincarities. We show that strongly 
photoexcited graphene monolayers with 35fs pulses quasi-instantaneously build up 
a broadband, inverted Dirac fermion population. Optical gain emerges and directly 
manifests itself via a negative conductivity at the near-infrared region for the first 200fs, 
where stimulated emission completely compensates absorption loss in the graphene 
layer. To quantitatively investigate this transient, extremely dense photoexcited 
Dirac-fermion state, we construct a two-chemical-potcntial model, in addition to a 
time-dependent transient carrier temperature above lattice temperature, to describe 
the population inverted electronic state metastablc on the time scale of tens of 
femtoseconds generated by a strong exciting pulse. The calculated transient optical 
conductivity reveals a complete bleaching of absorption, which sets the saturation 
density during the pulse propagation. Particularly, the model calculation reproduces 
the negative optical conductivity at lower frequencies in the states close to saturation, 
corroborating the observed femtosecond stimulated emission and optical gain in the 
wide near-infrared window. 
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1. Introduction 

Despite the well-established linear optical properties in graphene, which is marked 
by a universal absorption A = ira = 2.3% ranging from near-infrared to visible 
light, [U El [3] significantly less attention has been paid to the ultrafast nonlinear 
optical properties. Important for future photonic and optoelectronic applications, [1] 
carrier dynamics in graphene after been driven far out of equilibrium needs to be 
understood. Still, ultrafast spectroscopy studies are reported recently to show unusual 
properties. [61 El El El HU1 HU [121 HH1 HH HS1 HS1 HZ1 HH [191 EQl [21] Especially, the 
observation of nonlinear absorption when applying an ultrashort intense laser pulse 
to monolayer graphene reveals an extremely dense, quasithermal photoexcited-carrier 
state created by strong pumping on 10 fs time scale and metastable for several tens 
of femtoseconds, which implies a unique transient electronic state in the important 
emerging material graphene. |18| 

When strongly driven out of equilibrium by coherent light, the excited carriers 
subsequently participate in several dynamical processes during relaxing back to its 
equilibrium. Among them are the carrier decoherence, thermalization, cooling, and 
electron-hole recombination. If the excitation pulse is short enough, by observing 
the responses followed from right after the pump, we can identify the typical time 
scales associated with theses processes. Facilitated by recent ultrafast spectroscopy 
measurements, certain progress has been made. It is recognized that the ultrafast carrier 
dynamics in graphene is different from that in metals and semiconductors. First of all, 
dimensional estimates for Dirac fermions in graphene yield the carrier decoherence time 
r dc ~ 7^y^.[22] which becomes rather short, r^ c ~ lfs, for pump photon energy on the 
order of leV. Rapid carrier-carrier scattering gives rise to an electronic thermalization 
time Tth on the order of 10 fs. [T2| ITS] Time-resolved studies observe different cooling time 
scales with the shortest r c ~ lOOfs from electron-optical phonon coupling and longer 
electron-hole recombination time r r on the ps time scale. [51 [HI [TOl [61 [121 [131 [231 12H [22] 
Thus, when applying an ultrashort pump pulse of r p ~ lOfs, the distinct time scales 
Tth T c r entails an intermediate electronic state purely determined by carrier-carrier 
scattering. Unlike in most semiconductors where r t h ~ lOOfs > r p , the extremely short 
decoherence time and the rapid carrier-carrier scattering quickly deplete the phase space 
from the neighbourhood of optical excitation to the whole band to fully thermalize the 
carriers as illustrated in Fig. [1] (b). The analysis of Coulomb interaction between the 
Dirac- fermionic excitations shows a slow population imbalance relaxation, [221 I2E] in 
particular at high pulse intensity due to the suppression of Auger processes for strong 
pulse excitation, [25] although there is no gap. Thus for r t h < t < r Cjr the number of 
photoexcited carriers in each branch of the Dirac cone decays rather slowly. Therefore, 
the fast decoherence and thermalization together with slow imbalance relaxation support 
a unique population-inverted electronic state adiabatically formed after the strong 
pump excitation, a quasi-stable "hourglass" state on the time scale of some tens of 
femtoseconds, as shown in Fig. H(c). 
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Figure 1. Schematic illustration of the formation of population inverted electronic 
state in the intermediate time Tth < t < r c r . (a) Photoexcited carriers generated 
by ~10fs pump pulse; (b) The leading scattering processes of photo-excited carriers 
taking place in several femtoseconds: e + e — > e + e, h + h — > h + h, e + h — > e + h, 
which quickly establish individual thermalization in electron and hole carriers sharing 
a common electronic temperature T e due to the electron-hole scattering events, (c) 
After the internal thermalization, the photoexcited carriers form a population inverted 
hourglass-like electronic state characterized by two chemical potentials and a common 
electron temperature. 

In this article, we focus on this unique transient electronic state and provide details 
on ultrafast optical pump-probe spectroscopy measurement on monolayer graphene 
to reveal significant ultrafast optical nonlinearities, including nonlinear absorption 
saturation and near-infrared stimulated emission. These properties arise from a 
broadband, inverted Dirac fermion population induced by 35fs pulse excitation. Optical 
gain emerges and directly manifests itself via a negative conductivity at the near- 
infrared region for the first 100s of fs, where stimulated emission completely compensates 
absorption loss in the graphene layer. To quantitatively investigate this transient, 
extremely dense photoexcited Dirac-fermion state, we construct a simple model of 
a quasi-thermalized distribution with one electron temperature but two distinct 
chemical potentials associated with the electron- and hole-band, respectively. The 
characteristic electron temperature We find this transient state associated with high 
electron temperature T e up to 3000-4000K, which causes a broadband distribution 
extending to high energy that naturally explains the observed blueshifted component 
in photoluminescence spectrum. [7J El [9] We further explore the phase space capacity 
and identify the maximal photoexcitation density restricted by phase space filling, and 
individual chemical potential of each band are calculated. To understand the observed 
nonlinear optical behavior and the measured large saturation density, we calculate the 
optical conductivity for this nonequilibrium electronic state. The results show that 
the available phase space cannot be completely filled but will be saturated at a lower 
photoexcitation density due to the balance between absorption and emission. This 
calculated saturation density is in excellent agreement with the experimental value. [T8] 
Most interestingly, our model reproduces the experimental results that the nearly- 
saturated states created by a high-frequency pump are unstable to a low-frequency 
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pulse through stimulated emission to bring the system to a lower-level metastable state 
(illustrated in Fig. [2]), resulting in an optical gain phenomenon within the first 100s of fs 
after photoexcitation. The excellent agreement between theory and experiment further 
corroborates that our simple model captures the feature of the transient state at early 
timescale (< lOOfs) in the high excitation regime. Meanwhile, the comparison with the 
measured optical conductivity disfavors the equal-chemical-potential model to describe 
the transient states in graphene at the fs time scale, an outstanding issue debated in 
the community. 




Figure 2. Schematic illustration of (a) the femtosecond stimulated emission and 
optical gain when applying a low-frequency pulse to a high-density population-inverted 
state created by an intense high-frequency pump and (b) the resulting metastable low- 
density state. 

2. Experimental Details 

2.1. Spectroscopy measurement 

The experimental setup is shown in Fig. |3j In our experiment, the Ti:Sapphire 
amplifier with center wavelength 800nm, pulse width 35fs at 1kHz repetition rate is 
used. This further drives an optical parametric amplifiers tunable with tunable optical 
pump pulses covering 572-2400 nm allowing for both degenerate and non-degenerate 
pump/probe spectroscopy. The laser is further split into pump and probe paths. The 
pump beam, chopped as half harmonic of the laser repetition rate, directly excites the 
sample. The reflection of the probe beam, together with reference, is fed into an auto- 
balance detector, and the individual beams as well as the difference between them are 
picked up by three boxcar integrators. During the measurement, the pump fluence from 
few /xJ/cm 2 to mJ/cm 2 level is finely controlled. This way we can record pump- induced 
differential reflectivity changes AR/R with ~40fs time resolution and signal-to- noise 
ratio down to 5-10~ 5 . Similar experimental setups and details are described elsewhere, 
e.g., see Refs. [271 EEJ ESI EH ED [32] . 
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Figure 3. Experimental Schematics. HWP: half wavcplate; PB: polarizing beam 
splitter, FL: focal length. 



2.2. Samples 

Graphene was prepared from the thermal evaporation of SiC [33] with substrates used 
in the current experiments 6H-SiC(0001) purchased from Cree, Inc.. The samples were 
graphitized in UHV (P ~ 1 x 10~ 10 torr) by direct current heating of the sample to 
~1300°C measured with an infrared pyrometer with reading of the pyrometer adjusted 
to take account of the graphite emissivity reported in the literature [31]. The sample 
was not pretreated in a H2 atmosphere within a furnace which is a common practice 
because it excludes the formation of multi-step heights and easier control of thickness. 
The layer thickness (whether single layer Gl or bilayer G2) was controlled by the heating 
rate: faster one-step heating rates (within 2-3 seconds to reach 1300°C) result in large 
Gl domains while multiple heating steps with a slower rate (30seconds to reach 1350°C) 
result in samples with large G2 areas. Graphene thickness was identified using contrast 
thickness [35] and with step heights changes between different regions which were found 
to be combinations of only two steps ,i-e., 0.25nm (of SiC), and 0.33nm (of graphene ) as 
explained in Ref.[33]. Fig. lc in Ref. [H] shows a large 2/imx 2/xm (left) Gl formed after 
heating with the fast rate. The atomic scale image is shown to the right with the lxl 
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unit cell seen with lattice constant 0.246nm and intensity modulation due to the 6\/3 is 
also seen. The tunneling conditions are -0.5V, InA. The high intensity of the modulation 
and the resolution of the 6 atoms of the graphene ring indicate that this is predominantly 
Gl(in excess of 90% of the area). The detailed growth conditions, characterizations and 
doping (~0.4 eV for monolayer) of the obtained epitaxial graphene on SiC are extensively 
established by our papers [331 ESj and many others in the literature, e.g., [371 Ell 139] . 



3. Threshold reflection coefficient and optical gain 

Considering graphene on a substrate with dielectric constant e s , the amplitude of the 

reflected and transmitted waves for a normal incident beam follow from Maxwell's 

equations along with the usual boundary conditions: 

1 — n s — Ana (cj) Ic 

r = 

1 + n s + Ana (u) jc 

2 

t = 1 

1 + n s + Ana (u) /c 

Here a (u) is the complex optical conductivity. The common reflection and transmission 
coefficients are determined by R = \r\ 2 and T = n s \t\ 2 . In case a = holds that 
R+T = 1. The presence of a finite conductivity in the graphene sheet leads to absorption 

A g = ±(l + n.) 2 (l-T-R). (2) 

where it is custom [30] to introduce the coefficient (1 + n s ) 2 /A such that A g corresponds 
to the absorption coefficient of a suspended graphene sheet. 

Following Ref.jlQ] we can introduce the reflection of the substrate (for a = 0) 



l + n t 



and of the substrate with graphene R, 



+9 



Rs+g 



1 — n s — Ana (cu) jc 



2 



(4) 



1 + n s + Ana (u) jc 

For the complex optical conductivity of graphene in equilibrium and at T = it holds 

^> r =°> = s(^-^-i ln (^l) 2 )- (5) 

Near the jump in the optical conductivity at u = 2yU, the imaginary part of the 
conductivity has a logarithmic divergence which is smeared out in case of finite 
temperatures. Since a (u) is of order e 2 /h~, it holds that a (u) jc is of order of the 
finestructure constant of quantum electrodynamics oqed = e 2 / (he) w 1/137 1. This 
allows for an expansion in a (u) / c. It follows 
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Thus, the reflection coefficient to leading order in ojqed is fully determined by the real 
part of the optical conductivity a' (u) = Recr(cu), the imaginary part only enters at 
higher orders. For the transmission and absorption coefficients follows in the same limit 

4n s 8n s 4vr n( 2 \ 

1 = k o — o (oj) + U a opn 

(l + n s f (n s + l) 3 c K J y 

Air 

A 9 = -a' (u) + O (a| ED ) . (7) 



c 

This yields the result 



of Ref. 



Eqj6] enables us to determine a threshold value for the reflectivity that corresponds 
to a negative optical conductivity and thus to a behavior with optical gain. From Eqj6] 
it follows for the reflection after delay time r: 

R s+g (r) - R s+g (0) 



AR/R = 



R, 
in 



/W-^'(O)) (9) 
nj — 1 c 

Using the experimentally established value a' (0) = e 2 / (4h) for the optical conductivity 
prior to the pulse, it follows that a' (r) < if AR/R < AR/R\ c where 

Afl/flL = -5^. (10) 

With n s = 2.7 it follows AR/R\ c = —1.4582%. If for some reason the dielectric constant 
of the substrate is larger that 2.7, this would only reduce the critical value of AR/R and 
we would only underestimate the regime where a < 0. Given that our data yield the 
magnitude of AR/R as big as 1.9%, it follows that we have o < as long as n s > 2.41. 
In the literature, the uncertainty of n s = 2.7 is ±0.1. The smallest index of SiC is 2.55 
in the THz range. These results demonstrate that our conclusion o < is robust. 

In Fig. HI we experimentally determine the existence and value of the threshold 
AR/R\ c = —1.4582% for zero conductivity in our sample. This further demonstrates 
unambiguously that the reflectivity geometry in current sample provides a direct 
measurement of the real part of conductivity a of the graphene layer (or absorption), 
which directly accessing the gain/loss processes. This also demonstrates again our 
sample is graphene monolayer, consistent with conclusion from STM. 
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Figure 4. The differential reflectivity determined by the measurements with R g + S 
or without R s graphene monolayer on SiC substrate. The reflection from the zero 
conductivity in pumped garphene/SiC exactly corresponds to the case of bare SiC 
substrate. Consequently, the threshold AR/R\ C for zero conductivity can be directly 
determined from the curve, which is consistent with value used in the manuscript 
—1.46%. 
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Using the same reasoning we can relate the reflectivity to the absorption coefficient 

n — 1 4 

R s+9 (r) = R s + -— 2 A (r) (11) 
n s + 1 (1 + 71 1 



and obtain 

A g (r) - A gfi _ R s+g (r) - R s+g (0) - 1 + 4A g , 



^9,0 (0) 4A 9i0 

that will be used in our analysis of the density of transient electrons and holes, where 
A gfl = A g (r = 0) = 7ra QE D- 
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4. Stimulated infrared emission and optical gain 
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Figure 5. (a) Ultrafast AR/R at 1.55 eV pump and 1.33 eV probe for different 
pump fuences from 480 to 4759/1 J/cm 2 (b) Experimental values (rectangles) vs. theory 
(line) for peak transient conductivity, showing negative conductivity above a threshold 
pump flucncc. (c) The extracted peak transient transmission as function of the pump 
fluence clearly shows the positive transmission change, nonlinear saturation, and that 
the critical vaule ^?-| c (blue line) for zero conductivity indeed occurs. 

Here we provide a set of pump fluence dependence data at probe photon probe energy 
at 1.33 eV, as shown in Figs. [5^l and [5b. Our conclusions w.r.t. stimulated emission and 
optical gain are based on the observed negative conductivity in strongly photoexcited 
graphene, which is fully consistent with the complementary data presented in Ref [18J. 
In addition, following the similar analysis Eqs. (6) and (7), the differential transmission 
of our sample can be extracted by the information of the differential reflectivity or the 
subsequently derived conductivity of the graphene sample. More importantly, there 
also exists a threshold value for the photoinduced differential transmission -f-\ c that 
corresponds to a zero optical conductivity, above which the optical gain has to emerge 
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because of the negative conductivity. 

AT 2na.QED 



(13) 



T n s + l 

From Eqj6] and EqJT] it follows for the differential transmission after delay time r: 

AT AR 1 — n s 

-W = -g-M ■ — . (i4) 

The extracted peak transient transmission as function of the pump Fluence, as shown 
in Fig. [5b, clearly shows the positive transmission change and, mostly critically, that 
the critical value (blue line) for zero conductivity indeed occurs. While those data for 
the transmission were obtained indirectly, from our reflectivity measurements, a direct 
measurement of the transmission would be an important confirmation of our results 
ideally using large area free standing graphene monolayer samples. 

5. Analysis of the density of transient electrons/holes 

The amplitude of the time dependent absorption A as function of pump fluence can be 
derived from the measured differential reflectivity by applying the Fresnel equations in 
thin film limit @Dl ED] 

AACg = AR g+8 (I p ) ng-l + 4A g , 

Ao Rg+s 4A 9j0 

where Rg+s and AR g+s are the static reflectivity and pump-induced reflectivity changes 
for the graphene monolayer (g) on the substrate (s) with index n s = 2.7. A 9y o is 
the absorption of graphene monolayer without pump, which takes a universal value 

2 2 

of Agfi = nj- ~ 0.023, as determined by the universal a.c. conductivity <r = fr. 
Here the AA(I p )/Aq is the relative differential absorption of graphene on the substrate: 
A) = t^wA^o, which yields ^ = ^a. Therefore Eq.flED follows from Eq.(Q2}. 
The peak amplitude A(I P ) = A + AA(I p ) gradually diminishes as increasing the pump 
fluence. From the measured transient saturation curve above, one can extract the density 
of photoexcited electrons(holes) in graphene after the propagation of a single laser pulse 
of 35 fs (t p ) with pump fluence I p 

r°° 1 r°° dt 

n cx (I P )= -n ex (t,I p ) = — —I(t,I p )A(t), (16) 

J-oo T P J -oo T p 



where I(t, I p ) is the Gaussian pulse envelop I(t, I p ) = I p y exp 
such that the total pulse fluence is I p = J^^I^t, I p ). Since A(t) = A + AA(t) 



-4 In 2 ^2 

V 



normalized 



A (l + ^f), we have 



I(t,I p )A ( AA(t)\ 
n cx (t,/ p ) = — ^_M+_j. (17) 
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Applied to graphene where r th <C r p , A (t) is determined by the adiabatic dependence of 
the absorption on the pump fluence with /partial {t,Ip) = f_ y-I(tf,I p ). Consequently, 
EqJTTl becomes 



I(t,I p )A ( AA(I partial (t,I p )) 



(18) 



hco \ A 
We determine A (/partial) experimentally from the reflectivity data of Fig. [6^, combined 
with Eq]To"j as discussed above. Finally, from Eqs. (I15p - (jl8p we have 



^■ex (/p) 



I(t, I P )A 



x 



ARg+s (/partial (ty Ip)) 

Ro 



nz 



4A 



9,0 



4A 



dt 



(19) 



The result is shown in Fig. Eb, which clearly shows that using the actual absorption 
A (£), instead of A , is crutial to understand the high density regime of fs dynamics in 
graphene discussed here. 
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Figure 6. (a) The peak ^ : \ pea k as function of the pump fluence (black squares) 
measured by degenerate differential reflectivity at 1.55 cV for the graphene monolayer 
and the corresponding conductivity change (red solid dots). Blue arrow marks the 
threshold for zero conductivity AR/R\ C = —1.4582% (see text). Dashed line: linear 
dependence (guide to the eyes), (b) The extracted transient fermion density at 40 fs 
(blue dots), as explained in the text, which is significantly lower than Aolp/huj obtained 
from the universal conductivity (open circles), as illustrated in shadow area. 



6. Theory 

6.1. Model 

We construct the theoretical model based on the following considerations: the above 
analysis of the time scales associated with different dynamical processes shows that the 
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internal thermalization time r t h of the electron system is much shorter than the cooling 
r c and the electron-hole recombination time r r ; Simultaneous conservation of momentum 
and energy for Dirac-fermion scattering indicates the leading scattering processes are 
e + e — > e + e, h + h — > h+ h, e + h-^-e + hm strong excitation regime with suppression 
of Auger processes. [25] For the intermediate time regime r t h < t < r c , T , i-e., after 
absorbing photons but before losing energy to the lattice, photoexcited electrons and 
holes in graphene quickly establish separate thermal equilibrium through carrier-carrier 
scattering. This gives rise to sharply separated chemical potentials for the two bands, 
similar to the case of graphite thin film. [12] Due to the scattering between electron- and 
hole-carriers, e + h —> e+h, the whole electronic system obtains a common temperature, 
the electron temperature that differs from the lattice temperature. Therefore we 
characterize this intermediate electronic state by two Fermi-Dirac distributions 

/±(*) = — i (2°) 



exp 



£ ±( fc )~A f ± 
k B T e 



for the upper (+) and lower (— ) bands of the Dirac spectrum e±(k) = ±hvk with 
separate chemical potential fi± but the same electron temperature T e . Note that for 
the convenience of theoretical calculation, we use the upper- and lower-band electron 
picture instead of the electron- hole picture, which are related via /+ = f e , /_ = 1 — f\, 
and yU + = fi e , \i- = —fih- To solve for electron temperature and chemical potentials, 
we take into account the following conservation laws: 1) the total number of electrons 
before and after pump excitation is the same, 2) a pseudo-conservation law due to the 
slow population imbalance relaxation valid in strong excitation regime: the photoexcited 
carrier number in the intermediate state stays the same as that of right after the pump 
excitation, 3) and the above described adiabatic process requires that the absorbed 
photon energy is kept in the electron system until the formation of the quasi-thermal 
distribution ( |20|) . These three conditions are expressed as 



n to t = n + + n_ = n°, + n° 



??,, 



n+ 



n. 



n 



71- 



(21) 
(22) 



n ex hu = u — u , (23) 

where n tot = N tot /L 2 represents the total density of electrons in the system, n cx = 
N ex /L 2 refers to the density of photoexcited carriers, n± {n±) indicate the electron 
densities in the intermediate (initially equilibrium) state, and u (u°) represents the 
intermediate (initial) energy density of the whole electron system while hu is the pump 
photon energy. Applying the distribution ( 1201) to Eqs. (l2~TT) -(l23l) and taking into account 
the valley and spin degeneracy in graphene, we obtain the following expressions in terms 
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**± 



z± = e k B T e with initial temperature T = 300K, 



S 



9 (fc B T e ) 2 
2tt (/w) 2 
9 (k B Ty 

2tt (n«) 2 



-Lio 



-Li 2 (- 



-2;+ + 



Li-; 



Li 9 - 



(24) 



with 5 referring to the initial doping density with respect to the neutrality point, and 
9 1 



2tt (hv) 2 

9 1 

2vr (fa;) 2 



{(fc B T e ) 2 [-Li 2 (-^ + )] 



-Li, 



(A: B T ) 2 [-Li 
- (k B T ) 2 



*°)]} 



-Li, 



(25) 



as well as 



n ex hu 



9 (k B T e 

TV (hv)' 

9 (k B T 
(hv) 



-Li 3 (-z+) - Li 3 



0\3 



7T 



-^3 H°) - Li 3 (~) 



(26) 



where g = 4 is the flavor index taking into account the valley and spin degeneracies, v 
represents the Fermi velocity v ~ 1.1 x 10 6 m/s (355c) in graphene, and the polylogarithm 
is defined by a power series IA s (z) = X^Li S"- Solving the three equations gives the 
transient electron temperature T e and the individual chemical potentials \i± = k B T e In z± 
at a given photoexcitation density n ex with initial temperature T° and initial chemical 
potential fi° associated with the equilibrium state before being excited. 

To perform numerical calculation, we introduce the dimensionless variables 

, n ex 8 k B T e fi ± hu 

/ex = — , x = -, t= — — , a±=—, U = (27) 
n n D D D 

with a choice for the upper momentum cutoff A to define the energy scale D = hvA 
and the density scale n = — . Here we choose A such that 7rA 2 = |(27r) 2 /y4 where 
A = 3 3/2 ag/2 is the area of the hexagonal unit cell. Note that these dimensionless 
units are solely introduced for computational convenience. None of our final expressions 
depends on the actual values of A, D or n, as these quantities cancel in the final results 
(see for example Eq. (1341) ). In terms of the dimensionless variables, the equations are 
expressed as 



tot 



1 



X 



—Lie 



Lig 



9 ma 2 



-Li 2 (-z°) + Li 2 



(28) 



and 



/ex = f{t 2 [-Li 2 (-, + )]-(t°) 2 [-Li 2 (-Z )]} 



-Li 2 ( 



o\2 



(*°) 



" Li 2 I 'To 



(29) 
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as well as 



-Li 3 (-z+) - Li 3 



0\3 



-(*•) 



-Li s (---") - Li 3 



In the following analysis, either equation set ( )24|) -( l26l) or the set (J2H 
employed for convenience. 



(30) 

M will be 



6.2. Characteristics of the Intermediate Electronic State 

We carry out our analysis for two cases: undoped graphene, i.e. the system at 
the charge neutrality point, and graphene on SiC substrate with a finite electron 
doping. It is easy to show that the hole-doped system is symmetric to the electron- 
doped system. By solving for the electron temperature T e and the chemical potentials 
/x± at different photoexcitation densities, we demonstrate the characteristics of this 
intermediate electronic state. 



6.2.1. Neutral System The simplest case is the system at the neutrality point, i.e. 
5 = 0, possessing particle-hole symmetry. Equation ( 1281) yields z° = 1 (fi° = 0) and 
z + = — = z (/i+ = — /!_ = fi), i.e., the lower-band chemical potential is always the 
opposite of the upper-band one in the neutral system. From Eq. (|29|) and ( 130]) we 
obtain the expression for the dimensionless temperature t 



and the relation 



> + (*°) 2 ff 
-Li 2 (-z) 



1/2 

(31) 



3/2 



with 



*«- V 4 J + OT 3 !C(3) (32) 

h(z) S tMz^l. (33 ) 
-Li 3 (-z) 

Since h(z) is monotonously increasing with an upper bound 3/y/2 in large z limit, it 
implies a maximum value of / ex : /<^ ax = ( 3 j-) ■ That is to say, there exists a phase 
space limit on the photoexcited carrier number: 



< ax = T^—y- (34) 



9 

167TV 2 

For instance, at hw = 1.55eV the phase space capacity becomes ^™ ax = 9.7779 x 

10 13 cm -2 . In this limit, the electron temperature approaches zero as z — > oo, 
t y (/,T X ) 1/2 = _3n , Q 

In 2 4 In z 
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For a 800nm pump with photon energy fku = 1.55eV, solving Eq. ( 13 ip and 
we obtain T e function of ri cx plotted in Fig. [7J 



T.(K) 



1000 




Figure 7. Plot of the electron temperature T e varying with photoexcitation density 
n ex in the neutral system. We can see that the system is rapidly heated up at lower 
densities, but is slowly cooled at higher densities. 



It shows that the electron temperature rises rapidly at low photoexcitation densities, 

but, instead of keeping heated up, it starts slowly dropping at higher densities and 

eventually approaches zero at a maximal density. 

The value of fi with respect to n ex is plotted in Fig. [HJ 



low density regime 



n ci (10 ,: cm" 2 ) 




n ex (10'W 2 ) 



Figure 8. Plot of the upper-band chemical potential /i + (n cx ), in units of cV, in the 
neutral system. The lower-band chemical potential /i_ = —fi+. Clearly, the upper- 
band chemical potential turns into negative at low densities before rising to positive 
at high densities. 



We clearly see that it turns into negative at low photoexcitation densities during 
the rapidly heating up time, but back to positive when temperature slowly decreasing. 

The down-turn behavior in electron temperature and the negative-to-positive 
transition in chemical potential signifies a crossover behavior that at small pump fluence 
the excited carriers form a hot and dilute classical gas, but with more carriers excited 
they gradually build up a quantum degenerate fermion system with temperature cooling 
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down in order to accommodate more electrons in the finite phase space. If the phase 
space could really be exhausted, the electron and the hole carriers would be pumped 
into zero temperature Fermi-Dirac distributions in which the carriers are closely packed 
with a sharp Fermi edge. 



6.2.2. Exhaustion of Phase Space Inspired by the analysis of the neutral system, we 
see that phase space capacity is exhausted at zero electron temperature. To obtain 
an analytical estimate of the maximal available phase space at different electron-doping 
levels, we assume the initial temperature to be zero for convenience. Equations f T2TT) - (T23T> 
are then simplified as 



5 



7T 



{hv) 



1 



0\2 



7T 



(hv) 



(35) 



7T 



(ho) 



'hw 



0\2 



1 



(/i™ ax )' 



7T (hvf 

3 / 0\3 



(/i+ ax ) + (-/^ max ) - 0°) 



(36) 



(37) 



37T (hv) 

In terms of the dimensionless variables defined in ( 127)) . we find a relation between the 
maximal photoexcitation density and the doping level from the above equations 



pax / f 



max 
cx 



X 



,3/2 



(./; 



max\3/2 _ 3/2 
cx I x 



(3? 



Solving this equation yields n™ ax at different doping densities as shown in Fig. [91 



n ex (max) 10 
(lO'W 2 ) 




5 10 15 

6 (io'W 2 ) 



Figure 9. Plot of the maximum photocxcited carrier density n™ ax , when the phase 
space is completely filled, as a function of initial doping density 5. The maximum 
photoexcitation density decreases with increasing doping level, as expected. 



As seen from formula (134)) for the neutral system, n™ ax ~ ^j, the available phase 
space rises with increasing pump frequency. On the other hand, phase space capacity 
decreases with increasing initial electron doping density, as expected. 
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However, for photoexcitation density Tl ex clS clll input parameter in our calculation, 
a critical question to ask is: does this maximal density n™ ax equal the saturation density 
n s e ^ in real pumping process, or in another word, can phase space be completely filled? 
We will answer this question later. 



6.2.3. Doped System Next we discuss the system away from the Dirac point with a 
finite electron doping, i.e., 5>0orx = |>0, which is often the case, e.g., in epitaxial 
graphene on SiC substrate. In this case, from Eqs. fl28l) - fl30l) we obtain the expression 
of t as a function of z + 

'-( » + S:TT - 

then the coupled equations are reduced to 

-LiJ-—) =C 1 [-U 2 (-z + )], (40) 



-Li 3 (-z+) - Li 3 



1 

Z- 



2/3 

= C 2 [-Li 3 (-2+)] , (41) 



with 

Ci(/ex) = 1 r > (42) 

U } / cx + 2(tO) 2 [-Li 2 (-^)]' 

_ (^f + (t°) 3 [-Li 3 (-,°)-Li 3 (-^]) 2/3 

C2(/cx) " if + (to) 2 [-Li 2 (-*o)] ■ (43) 

Solving the two equations (HP]) and (j4"Tl) we obtain z + and z_, which in turn gives t 
via Eq. ( l39l) . Finally, the physical quantities are derived through k^T e = tD, fi + = 
tD In z + , /!_ = t.D In z_ . 

To show the numerical results, we choose the experimental system of graphene on 
SiC substrate with an initial electron doping 6 = 1.17 x 10 13 cm~ 2 , corresponding to an 
initial chemical potential fi° = 0.4eV, being excited by the pump energy hu = 1.55eV. 
The phase space capacity is calculated from Eq. ( 1381) to be n™ ax = 8.34 x 10 13 cm -2 
at this doping level. The electron temperature T e (in units of Kelvin) changing with 
n ex (in units of 10 13 cm~ 2 ) is plotted in Fig.[TUJ Compared to the undoped system, 
the evolution of electron temperature with photoexcitation density is smoother and the 
electron temperature is lower due to the finite initial carrier density. 
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T.(K) 




2 4 6 8 



n ex (10 13 cm 2 ) 

Figure 10. Plot of the electron temperature T e varying with photocxcitation density 
n ex in epitaxial graphene on SiC substrate with initial chemical potential /i = 0.4eV. 
The non-monotonous behavior remains although the change of temperature with 
density is smoother compared to the neutral system. 



Figure [TT] shows the upper- and lower-band chemical potential fi + and /x_ (in units 
of eV) varying with n ex (in units of 10 13 cm~ 2 ). 



/U(eV) 




Figure 11. Plot of the upper- and lower-band chemical potentials /i + and in units 
of eV, varying with photocxcitation density n ex in epitaxial graphene on SiC substrate 
with initial chemical potential fx = 0.4cV. The upper-band chemical potential remains 
positive although it drops a little bit at low densities, while the lower-band chemical 
potential drops rapidly at rather low densities followed by increasing separation from 
the upper-band chemical potential. Clearly, the upper- and lower-band chemical 
potentials are not symmetric as in the neutral system. 



Clearly, due to the large initial electron-doping, the low-density classical gas phase 
for the upper-band electrons is absent now, although there is a tiny presence for the 
lower band. And in the doped case, the upper- and lower-band chemical potentials are 
not symmetric, /i + ^ —ft-, as they are in the neutral system. The separation between 
the two chemical potentials increases with photoexcitation density. 

6.2.4- Broadband Distribution and Blueshifted Photoluminescence A direct conse- 
quence to the high electron temperature T e ~ 3000 — 4000K and the slow population 
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relaxation is a broadband distribution of electron and hole excitations. This can be 
illustrated in the occupation number N e>h (e) = Z?(e)/ Cj / l (e), where D(e) = ,^\ 2 is the 
density of state at energy e and f e ,h{^) ar e the electron and hole distribution functions 
with f e = /+, fh = 1 — /_. Figure [T2l shows the electron and hole distribution at dif- 
ferent photoexcitation density in the neutral system and in the electron-doped system 
0° = 0.4eV). 



neutral system 

occupation 




-3 -2 -1 1 2 3 

energy (eV) 



doped system 




-3-2-10 1 2 3 

energy (eV) 

Figure 12. Electron (blue line) and hole (green line) occupation number, in units 
of n /ft v \2 1 distributed with respect to energy at different photoexcitation density 
states. We can see the broadband distribution of photoexcited carriers with a high- 
temperature tail up to almost 2 eV. 

The high temperature tail in the distribution extends the excited carriers well 
above the excitation energy 1.55eV up to 2-3eV. This coverage of higher energy states 
enables emission of photons with higher frequencies than that of the excitation photons, 
which exhibits blue-shift phenomena in the photoluminescence spectrum. This unusual 
blueshifted components have been observed in recent experiments [HEl Ej and is naturally 
explained in our model. 
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6.3. Optical Conductivity 

Now we are back to the question raised before. In order to identify the saturation density, 
we need to study the optical responses of this electronic state at different photoexcitation 
densities. In this section, we calculate the optical conductivity for the nonequilibrium 
intermediate state using Keldysh technique. By analyzing its behavior at high densities, 
interesting optical properties are revealed. 



6.3.1. General Formalism The low-energy noninteracting Hamiltonian of graphene can 
be written in the band representation as 



H ° = V ib [ E A HU(k)7a,A(k) 
a=l Jk X=± 



(44) 



with A = + or — corresponding to the upper or lower band and 7 ( [ A (k), 7 ai A(k) are 
the operators that create or annihilate a quasiparticle of flavor a (spin and valley) at 
the 2- dimensional wavevector k = (k x , k y ) in band A. In the band representation the 
current vertex becomes 

: /k kxe 2 \ 

Jk = ev I -a z — a y 1 , (45) 

where Pauli matrices due to the chiral structure of the Dirac fermions in 

graphene. Note that we set H = 1 during the derivation, but will recover it in the 
final results. 

In Keldysh formalism, the bubble diagram contributing to the optical conductivity 
gives the real part as 



Rea a a(uj) = — 



/doj d k [~ - - 

^2 Tr ja^A^uj' + uj)j^N v {uj') - ja^A^u^j^N^uj' + u) 

(46) 

for the direction a(/3) = x, y in the 2D graphene layer with the definitions 



4M = ^(^M-GrH), (47) 

-^<( W ). (48) 
Here the retarded and advanced Green's functions are matrices in band representation 

\ u±iO++ft+vk J 

where + (— ) sign associates with the retarded (advanced) Green's function, and the 
lesser Green's function is written as 
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2nif(e±(k) - li±)5(uj - e±(k) + /x), 



with fj, — h (fx + + ix J) and the Fermi function f(x) 



(51) 

2 . f-j — — — — jvj - e ./(fc B 1 T e)+1 - Note that distinct 
chemical potentials are employed in the distribution functions to characterize the 
nonequilibrium state but an average chemical potential is used in the spectral functions 
to avoid an artificial modification of the spectrum. 

Previous analysis (see Section 3) shows that the optical properties of graphene 
on an insulating substrate, such as reflection, transmission, and absorption, are fully 
determined by the real part of the optical conductivity Rea(uj) to leading order in 
the fine-structure constant of quantum electrodynamics «qed = §3 ~ "C 1. The 
imaginary part only enters at higher orders. Therefore, we make direct connections of 
the real part conductivity to the optical responses observable in experiments. 

From Eq.( l46j) we calculate the longitudinal conductivity which contains intraband 
and interband transitions. To show the transition processes specifically, introduce 
a>k,\(u) = 5(cu — Xvk + ft), and f\{uj) = f(co — XSfx). It follows the intraband and 
interband conductivities 



gn (ev) 



Real n ' ei (u 



gir (ev) 



dco'd 2 k 



(27T) 



COS' 



d^^M^ ) fl k,A(w +U) — (52) 



A=± 



UJ 



du'drk 



(27T) 



sin 



2_^a K \(uj )a k> x(o; +w) — (-5,3) 



A=± 



UJ 



where cos 9 



f, sin0=£. 



6.3.2. Intraband Transition First let us evaluate intraband conductivity, 
straightforwardly obtained from Eq. (??) 



It is 



Rea 



intra. 



9(ev) 
(2) 2 



-5{u) 



kdk 



dm 



dm 



In [(l + z + )(l + 



<9w 

k B T e S(uj) 



du 



-vk—fi- 



(54) 



where we have recovered the factor H on the last line. The delta-function will be 
replaced by a Lorentzian 6(u) — > aj2 _ l T( T -i)2 f° r further discussions, which is not our 
concern here. In equilibrium state z + 



becomes Reai 



Z- = z 
k B T° ln(2 + + e 



e^ , the intraband conductivity 
5(uj) which recovers the well- 

= 21n2^k B T°5{uj), 



intra , 



known expression in the neutral system at equilibrium Rea 
as expected. 

The most interesting observation from the intraband transition for the transient 
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electronic state is the modified Drude spectral weight 

e 2 

D = - In [(1 + z+){l + zZ 1 )] k B T e , (55) 

which can be significantly enhanced by the high electron temperature. But it could also 
be reduced when the reduction of chemical potential dominates at low densities in the 
neutral system. Here we show the Drude weight change, normalized by the equilibrium 
value D = \ In [(1 + z°)(l + 1/z )} k B T°, at different n ex in the neutral and electron- 
doped systems in Fig. [T3j 

The results for the neutral system exhibits a drop in Drude spectral weight at low 
densities due to the large drop in chemical potential as shown in Fig. [8] but quickly 
followed by large enhancement at higher densities. In the electron-doped system, Drude 
weight is always increasing but with much less enhancement than in the neutral system. 




12 3 4 

n„(10"cm : ) 



Figure 13. Drude spectral weigh D, normalized by the equilibrium value Dq, changes 
with photoexcitation density in neutral and doped (with initial chemical potential 
/i° = 0.4eV) system. In neutral system it drops at low densities due to the large 
drop in chemical potential as shown in Fig.|8l but followed by a large enhancement at 
higher densities. In the doped system, it is always enhanced but with a much smaller 
enhancement factor than that of the neutral system. 



6.3.3. Interband Transition In order to understand the optical response at high 
frequencies, as optical conductivity is dominated by interband transition for frequencies 
on the order of leV, we evaluate the interband conductivity from Eq.(??), 

Recxl ntcr/ 



9{e) 



(2) 4 
4H2 



fir 



2 



tanh 



hui + 2/i_ 
4k B T e 



tanh 



hui — 2/i H 



(56) 
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where the probe photon frequency u) > and we have reinserted the factor H on 



the last line. In equilibrium state, //_ 



Rea 



inter , 



4/12 



tann ^ 4fcsr0 



+ tanh 



neutral system in equilibrium Rea 



inter 



4k B T° 
~- fi tanh 



the interband transition becomes 
which gives the expression for the 



4fc s T° 



as expected. 



By studying the optical response to different photon energies at various 
photoexcitation densities, unusual optical properties of the transient electronic state 
are found, which will be discussed in the following. 



Femtosecond Absorption Saturation and Perfect Transparency Let us first consider the 
optical response to the pump frequency. An interesting observation from the interband 
transition formula (|56|) arises due to the two density-dependent chemical potentials. As 
shown in FigJH] and HU with increasing photoexcitation density, the separation of the 
two chemical potentials also gets larger. Then at a certain photoexcitation density such 
that 

fi + — //_ = huj, (57) 

the optical conductivity vanishes and the system becomes perfect transparent. For 
higher densities, fi+ — \i- > frou such that optical conductivity turns into negative, which 
implies a stimulated emission to keep the photoexcitation density from rising. This 
indicates that the absorption reaches zero and the number of excited carriers saturates 
at this density n s c ^ , which is stabilized by stimulated emission. 

To show the variation of optical conductivity with photoexcitation density, we again 
calculate for the neutral system and the electron-doped system with initial doping 
fjP = 0.4eV by applying pump photon energy at hu = 1.55eV. We show the results 
in Fig fl4l where they are normalized by the equilibrium value o"o = fr. 
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Figure 14. Intcrband conductivity a(hco = 1.55cV), normalized by the equilibrium 

2 

value (To = fw) varies with photoexcitation density n cx in the neutral and doped 
systems. It reaches zero, which indicates complete bleaching of absorption, at 
5.7 x 10 13 cm -2 in the neutral system and 4.8 x 10 13 cm~ 2 in the doped system, 
corresponding to the saturation densities. 

In both cases, the conductivity monotonously decreases due to the increasing 
electron temperature and separation of chemical potentials. The neutral system 
saturates at n^f (theory) = 5.7 x 10 13 cm~ 2 while the electron-doped system saturates at 
roughly n^* (theory) = 4.8 x 10 13 cm -2 . On the other hand, experimental measurement of 
the electron-doped system gives n|^(expt.) = 5.0 x 10 13 cm -2 , as shown in Fig. 6b, which 
is in excellent agreement with the theoretical value 4.8 x 10 13 cm -2 . This corroborates 
the correct description of the transient electronic state in our model. And it also answers 
the early posted question: the system saturates at a lower photoexcitation density before 
completely filling the available phase space, i.e., n s e ^ < n™ ax . 
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Figure 15. The calculated intcrband conductivity ct(oj), normalized by the equilibrium 
value (To = jf-, at higher (huj = 1.7eV, green line), the pump (hui = 1.55eV, magenta 
line), and lower (Huj = 1.2cV, blue line) frequency. It shows that in the vicinity of 
saturation density the optical conductivity at lower frequencies (but still high enough 
to mainly detect interband transition) becomes negative. This indicates a stimulated 
emission that enables an optical gain for the low- frequency probes. 

This theoretical calculation of an optical gain can serve as a test of our model, 
which have been used to simulate the optical differential reflectivity data performed in 
Fig. 5b (red line). The agreement is excellent. 

Femtosecond Stimulated Emission and Optical Gain It is easy to see that higher- 
frequency pump will saturate at higher density since it can open up more phase space. 
Then when applying a probe with frequency higher than the pump frequency, we will 
expect it can not detect the zero absorption, as long as it is still within the low-energy 
Dirac spectrum, as shown in Fig{15] (green line). However, if one applies a lower- 
frequency probe, but not too low such that it is still mainly detecting the change in 
interband transition, one would expect an optical gain in the vicinity of the saturation 
density, as the optical conductivity becomes negative in this regime as shown in FiglTol 
(blue line, yellow region). The appearance of negative conductivity signifies a stimulated 
emission that drives the system to a lower density as illustrated in FigfSJ We stress that 
the transient conductivity is negative in a regime below the pump frequency but above 
a certain frequency below which intraband transition becomes dominant. 
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Comparison of Two Model Calculations with Pump-probe Spectroscopy Measurements 
Here we further compare the calculated optical conductivity from the distinct-// model 
(/!+ 7^ /!_) discussed above and the equal-// model (//+ = //_) with the experimental 
value measured at probe photon energy hu> = 1.55eV and 1.16eV for a pump energy 
at fruj = 1.55eV. As shown in Figfl6l we compare the experimentally-extracted, 
transient conductivity at 40 fs [18] with the calculated conductivity as a function 
of the photoexcited carrier density n cx for two probe photon energies 1.55 eV and 1.16 
eV. The Fermi energy of the sample is ~ 0.4 eV. The model calculation with the 
distinct chemical potentials reproduces the salient features of the experiment including 
nonlinear saturation and optical gain. Excellent agreement between experiment and 
theory also demonstrates a faithful representation of the transient state at 40 fs by 
the model described in the manuscript. The model calculation with the same chemical 
potential clearly fails to account for the experimental observations. For the degenerate 
scheme, our theory (black dashed line) yields a — > and thus perfect transparency 
at n ex = 4.8 x 10 13 cm -2 . Once the system is driven into this regime, a balance 
between stimulated emission and absorption will lead to a transparency. For non- 
degenerate scheme by probing at 1.16 eV, our theory (black solid line) predicts a 
critical density 3.2 x 10 13 cm -2 for the transition from loss to gain. All of these results 
agree quantitatively with the experimental values 5.0 x 10 13 cm -2 and 3.4 x 10 13 cm -2 , 
respectively. 




Figure 16. Comparison of the calculation from the distinct- // model (// + 7^ M-) 
and the equal-// model (//+ = //_) with the experimentally measured transient optical 
conductivity at 40 fs after the 1.55eV pump at varying pump fiuence. The probe photon 
energies are Huj = 1.55eV (blue solid square) and 1.16eV (red solid square). Clearly, 
the distinct-// model calculation of the conductivity at 1.55eV (black dashed line) and 
1.16eV (black solid line) agrees quantitatively with the experimental data, in sharp 
contrast to the equal-// model results (green dashed line for 1.55eV and green solid line 
for 1.16eV). For the probe photon at lower frequency 1.16eV in the non-degenerate 
scheme, the transient conductivity becomes negative above a critical density exhibiting 
a transition from optical loss to gain, predicted by distinct-// model and substantiated 
by experiment. 
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7. Conclusions 

We have studied the electronic state in photoexcited graphene formed via rapid carrier- 
carrier scattering after strong photoexcitation but before energy relaxation that takes 
place on a longer time scale. We have provided evidence for the existence of pronounced 
femtosecond population inversion and broadband gain in strongly photoexcited graphene 
monolayers. These results clearly reveal the transient electron and hole potentials are 
separated on the time scale of 100s of fs. By characterizing the state in terms of two 
separate Fermi-Dirac distributions with a common electron temperature but distinct 
chemical potentials for the upper and lower bands, we showed that this intermediate 
electronic state is associated with high electron temperature T e up to 3000-4000K, which 
causes a broadband distribution extended to higher energy and is responsible for the 
observed blueshifted photoluminescence component. The analysis on the variation of 
electron temperature and chemical potentials with photoexcited carrier density in the 
neutral system clearly shows a crossover from hot dilute classical gas to dense quantum 
degenerate fermions. And unlike the phase space restriction in most semiconductors for 
a pump pulse on the order of lOfs, which is determined by the density of state at the 
optical excitation and the frequency width of the pulse, the fast depletion of phase space 
in graphene yields a broadband filling which significantly enlarge the accommodation of 
photoexcited carriers. 
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